Studying early structural changes in SOS1 mediated KRAS activation mechanism.

KRAS activation is known to be modulated by a guanine nucleotide exchange factor (GEF), namely, Son of Sevenless1 (SOS1). SOS1 facilitates the exchange of GDP to GTP thereby leading to activation of KRAS. The binding of GDP/GTP to KRAS at the REM/allosteric site of SOS1 regulates the activation of KRAS at CDC25/catalytic site by facilitating its exchange. Different aspects of the allosteric activation of KRAS through SOS1 are still being explored. To understand the SOS1 mediated activation of KRAS, molecular dynamics simulations for a total of nine SOS1 complexes (KRAS-SOS1-KRAS) were performed. These nine systems comprised different combinations of KRAS-bound nucleotides (GTP/GDP) at REM and CDC25 sites of SOS1. Various conformational and thermodynamic parameters were analyzed for these simulation systems. MMPBSA free energy analysis revealed that binding at CDC25 site of SOS1 was significantly low for GDP-bound KRAS as compared to that of GTP-bound KRAS. It was observed that presence of either GDP/GTP bound KRAS at the REM site of SOS1 affected the activation related changes in the KRAS present at CDC25 site. The conformational changes at the catalytic site of SOS1 resulting from GDP/GTP-bound KRAS at the allosteric changes may hint at KRAS activation through different pathways (slow/fast/rare). The allosteric effect on activation of KRAS at CDC25 site may be due to conformations adopted by switch-I, switch-II, beta2 regions of KRAS at REM site. The effect of structural rearrangements occurring at allosteric KRAS may have led to increased interactions between SOS1 and KRAS at both the sites. The SOS1 residues involved in these important interactions with KRAS at the REM site were R694, S732 and K735. Whereas the ones interacting with KRAS at CDC25 site were S807, W809 and K814. This may suggest the crucial role of these residues in guiding the allosteric activation of KRAS at CDC25 site. The conformational shifts observed in the switch-I, switch-II and alpha3 regions of KRAS at CDC25 site may be attributed to be a part of allosteric activation. The binding affinities, interacting residues and conformational dynamics may provide an insight into development of inhibitors targeting the SOS1 mediated KRAS activation.

MolToxPred: small molecule toxicity prediction using machine learning approach.

Different types of chemicals and products may exhibit various health risks when administered into the human body. For toxicity reasons, the number of new drugs entering the market through the conventional drug development process has been reduced over the years. However, with the advent of big data and artificial intelligence, machine learning techniques have emerged as a potential solution for predicting toxicity and ensuring efficient drug development and chemical safety. An ML model for toxicity prediction can reduce experimental costs and time while addressing ethical concerns by drastically reducing the need for animals and clinical trials. Herein, MolToxPred, an ML-based tool, has been developed using a stacked model approach to predict the potential toxicity of small molecules and metabolites. The stacked model consists of random forest, multi-layer perceptron, and LightGBM as base classifiers and Logistic Regression as the meta classifier. For training and validation purposes, a comprehensive set of toxic and non-toxic molecules is curated. Different structural and physicochemical-based features in the form of molecular descriptors and fingerprints were employed. MolToxPred utilizes a comprehensive feature selection process and optimizes its hyperparameters through Bayesian optimization with stratified 5-fold cross-validation. In the evaluation phase, MolToxPred achieved an AUROC of 87.76% on the test set and 88.84% on an external validation set. The McNemar test was used as the post-hoc test to determine if the stacked models' performance was significantly different compared to the base learners. The developed stacked model outperformed its base classifiers and an existing tool in the literature, reaffirming its better performance. The hypothesis is that the incorporation of a diverse set of data, the subsequent feature selection, and a stacked ensemble approach give MolToxPred the edge over other methods. In addition to this, an attempt has been made to identify structural alerts responsible for endpoints of the Tox21 data to determine the association of a molecule with a plausible downstream pathway of action. MolToxPred may be helpful for drug discovery and regulatory pipelines in pharmaceutical and other industries for in silico toxicity prediction of small molecule candidates.

Structural insights into the activation mechanism of phosphoinositide 3-kinase alpha.

Phosphoinositide 3-kinases (PI3Ks) are lipid kinases known to regulate important cellular functions by phosphorylating the inositol ring of inositol-phospholipids (PtdIns) at 3' position. The PI3Kα is a heterodimer and the activation of the catalytic subunit (p110α) is regulated by its regulatory subunit (p85α). The current work deals with studying the activation mechanism of the PI3Kα using multi micro-second molecular dynamic simulations. Structural changes involved in activation mechanism is studied by gradually releasing the inhibitory effects of different domains of regulatory subunit namely, n-terminal SH2 (nSH2) and inter SH2 (iSH2). The observation shows that even in the presence of n-terminal and inter SH2 domain (niSH2) of regulatory subunit, the catalytic domain has some intrinsic activation activity and the presence of c-terminal SH2 (cSH2) domain may be required for complete inhibition. The release of nSH2 domain leads to loss of interactions between iSH2 domain (regulatory subunit) and C2 and kinase domain (catalytic subunit). The study shows that early events in the activation mechanism involve the movement of the ABD domain of the catalytic subunit along with the linker region between ABD and RBD region which may lead to movement of ABD closer to the CLobe of the kinase domain. This movement is essentially as it triggers the rearrangement of CLobe especially the catalytic loop and activation loop which bring catalytic important residues closer to ATP and PIP2(phosphatidylinositol-4,5-bisphosphate). Water mediated interaction analysis reveal that water may be playing an important role in the transfer of phosphate from ATP to PIP2. The study shows that initial signal for release of inhibitory effect of the regulatory subunit might be propagated through the linker region between ABD and RBD through allosteric effect to different regions of the protein. These understanding of early events during the activation mechanism may help in the design of better therapeutic targeting PI3K.

Novel alkyl-substituted 4-methoxy benzaldehyde thiosemicarbazones: Multi-target directed ligands for the treatment of Alzheimer's disease.

Alzheimer's disease (AD) is a devastating neurodegenerative disorder affecting mental ability and interrupts neurocognitive functions. Treating multifactorial conditions of AD with a single-target-directed drug is highly difficult. Thus, a multi-target-directed ligand (MTDL) development strategy has been developed as a promising approach for the treatment of AD. Herein, we have synthesized two novel thiosemicarbazones as MTDLs and reported their bioactivities against diverse neuropathological events involved in AD. In vitro studies revealed that both compounds exhibited promising anticholinesterase activity (AChE, IC50 = 15.98 µM, MZET and IC50 = 30.23 µM, MZMT), well supported by a detailed computational study. Both analogs have shown good thermodynamic behaviour and stability through interactions with characteristic amino acid residues throughout simulation of 100 ns against acetylcholinesterase enzyme. In an electrophysiology assay, these analogs have shown a characteristic inhibitory response against the GluN1-1a + GluN2B subunit of N-methyl-D-aspartate receptors. Pre-treatment of BV-2 microglial cells with MZET effectively decreased nitrite production compared to nitrite produced by lipopolysaccharide-treated cells alone. Further, the effect of MZMT and MZET on autophagy regulation was determined using stably transfected SH-SY5Y neuroblastoma cells. MZET significantly enhanced the autophagy flux in neuroblastoma cells. A significant decrease in copper-catalysed oxidation of amyloid-ß in presence of synthesized thiosemicarbazones was also observed. Collectively, our findings indicated that these analogs have potential as effective anti-AD candidates and can be used as a prototype to develop more safer multi-targeted anti-AD drugs.

Understanding the binding affinities between SFRP1CRD, SFRP1Netrin, Wnt5B and frizzled receptors 2, 3 and 7 using MD simulations.

cWnt-signalling plays a crucial role in stem cell maintenance and tissue homeostasis. Secreted frizzled-related proteins(SFRP), Wnt inhibitors consist of the N-terminal cysteine rich domain(CRD) and the C-terminal netrin(NTR) domain. SFRP1 binds to the Wnt ligands and frizzled receptors(FZ) either through its SFRP1CRD or through its SFRP1Netrin domains; however, very little is known on these binding affinities. Here, we attempted to understand the interactions and binding affinities of SFRP1-Wnt5B, SFRP1-FZ(2, 3 & 7) and Wnt5B-FZ(2, 3 & 7) that are mainly expressed in murine hair follicle stem cells. SFRP1CRD, SFRP1Netrin, Wnt5B and FZ(2, 3 & 7) structures were built using homology modelling, followed by their molecular dynamics simulations. SFRP1CRD showed lower fluctuation when in complex with FZ2, FZ3 and FZ7 and Wnt5B as compared to SFRP1Netrin using RMSF and RMSD. However, free energy showed SFRP1Netrin was energetically more stable than SFRP1CRD. SFRP1Netrin formed more number of interactions with FZ as compared to SFRP1CRD. Importantly, SFRP1Netrin favoured binding to the FZ receptors(FZ3 > FZ7 > FZ2) as compared to Wnt5B ligand. Conversely, the SFRP1CRD showed more affinity towards the Wnt5B ligand as compared to FZ receptors. Wnt5B showed the best binding affinity with FZ3 followed by SFRP1CRD and SFRP1Netrin. Therefore, SFRP1Netrin can bind to the FZ3 with higher binding affinity and may inhibit non-canonical Wnt-signalling pathway. Our study provides the comprehensive information on the binding affinities among the Wnt5B, SFRP1CRD/Netrin and FZ(2, 3 & 7). Thus, this information might also help in designing novel strategies to inhibit aberrant Wnt-signalling.Communicated by Ramaswamy H. Sarma.

Structural insight into the binding interactions of NTPs and nucleotide analogues to RNA dependent RNA polymerase of SARS-CoV-2.

RNA dependent RNA polymerase (RdRP) from positive-stranded RNA viruses has always been a hot target for designing of new drugs. Major class of drugs that are targeted against RdRP are nucleotide analogues. Extensive docking and molecular dynamics study describing the binding of natural nucleotides (NTPs) and its analogues leading to significant structural variation in the RdRP has been presented here. RdRP simulations in its apo, NTP-bound, and analogue-bound form have been performed. Nucleotide analogues included in this study were, favipiravir, galidesivir, lamivudine, ribavirin, remdesivir and sofosbuvir. The conformational flexibility of the RdRP molecule has been explored using principal component (PCA) and Markov state modeling (MSM) analysis. PCA inferred the presence of correlated motions among the conserved motifs of RdRP. Inter-domain distances between the finger and thumb subdomain flanking the nascent RNA template entry site sampled open and closed conformations. The ligand and template binding motifs F and G showed negatively correlated motions. K551, R553, and R555, a part of motif F appear to form strong interactions with the ligand molecules. R836, a primer binding residue was observed to strongly bind to the analogues. MSM analysis helped to extract statistically distinct conformations explored by the RdRP. Ensemble docking of the ligands on the Markov states also suggested the involvement of the above residues in ligand interactions. Markov states obtained clearly demarcated the open/closed conformations of the template entry site. These observations on residues from the conserved motifs involved in binding to the ligands may provide an insight into designing new inhibitors.Communicated by Ramaswamy H. Sarma.

Molecular dynamics of hERG channel: insights into understanding the binding of small molecules for detuning cardiotoxicity.

Evaluation of cardiotoxicity potential of new chemical entities (NCEs) has lately become one of the stringent filters in the drug discovery and development process. Cardiotoxicity is caused mainly by the inhibition of human ether-a-go-go related gene (hERG) channel protein. Inhibition of the hERG channel leads to a life-threatening condition known as cardiac arrhythmia. Knowledge of the structural behaviour of the hERG would aid greatly in the design of new drug molecules that do not interact with the protein and add to the safety index. In this study, a computational model for the active-state of hERG was developed. This model was equilibrated by performing the molecular dynamics simulations for 100 ns followed by clustering and selection of a representative structure based on the largest populated cluster. To study the changes in the protein structure on inhibition, three inhibitory ligands, namely, dofetilide, cisapride and terfenadine were docked, followed by molecular dynamics simulations of 200 ns for the apo and each ligand-bound structure. It was observed that docking and simulation studies of the hERG model exhibited noticeable conformational changes in the protein upon ligand-binding. A significant change in the kink of the S6-transmembrane helix was observed. Inter-chain distances between the crucial residues Y652 and F656 (present below the ion-selectivity filter), their side-chain orientation and hydrogen bonding indicated a probable collapse of the pore. These changes may infer the initiation in transition of hERG from an open to an inactive state. Hence, these findings would help in designing compounds devoid of hERG inhibition with reduced cardiotoxicity.Communicated by Ramaswamy H. Sarma.

Natural plant products as potential inhibitors of RNA dependent RNA polymerase of Severe Acute Respiratory Syndrome Coronavirus-2.

Drug repurposing studies targeting inhibition of RNA dependent RNA polymerase (RdRP) of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) have exhibited the potential effect of small molecules. In the present work a detailed interaction study between the phytochemicals from Indian medicinal plants and the RdRP of SARS-CoV-2 has been performed. The top four phytochemicals obtained through molecular docking were, swertiapuniside, cordifolide A, sitoindoside IX, and amarogentin belonging to Swertia chirayita, Tinospora cordifolia and Withania somnifera. These ligands bound to the RdRP were further studied using molecular dynamics simulations. The principal component analysis of these systems showed significant conformational changes in the finger and thumb subdomain of the RdRP. Hydrogen bonding, salt-bridge and water mediated interactions supported by MM-GBSA free energy of binding revealed strong binding of cordifolide A and sitoindoside IX to RdRP. The ligand-interacting residues belonged to either of the seven conserved motifs of the RdRP. These residues were polar and charged amino acids, namely, ARG 553, ARG 555, ASP 618, ASP 760, ASP 761, GLU 811, and SER 814. The glycosidic moieties of the phytochemicals were observed to form favourable interactions with these residues. Hence, these phytochemicals may hold the potential to act as RdRP inhibitors owing to their stability in binding to the druggable site.

An insight into the inhibitory mechanism of phytochemicals and FDA-approved drugs on the ACE2-Spike complex of SARS-CoV-2 using computational methods.

The S-glycoprotein (Spike) of the SARS-CoV-2 forms a complex with the human transmembrane protein angiotensin-converting enzyme 2 (ACE2) during infection. It forms the first line of contact with the human cell. The FDA-approved drugs and phytochemicals from Indian medicinal plants were explored. Molecular docking and simulations of these molecules targeting the ACE2-Spike complex were performed. Rutin DAB10 and Swertiapuniside were obtained as the top-scored drugs as per the docking protocol. The MD simulations of ligand-free, Rutin DAB10-bound, and Swertiapuniside-bound ACE2-Spike complex revealed abrogation of the hydrogen bonding network between the two proteins. The principal component and dynamic cross-correlation analysis pointed out conformational changes in both the proteins unique to the ligand-bound systems. The interface residues, His34, and Lys353 from ACE2 and Arg403, and Tyr495 from the Spike protein formed significant strong interactions with the ligand molecules, inferring the inhibition of ACE2-Spike complex. Few novel interactions specific to Rutin-DAB10 and Swertiapuniside were also identified. The conformational flexibility of the drug-binding pocket was captured using the RMSD-based clustering of the ligand-free simulations. Ensemble docking was performed wherein the FDA-approved database and phytochemical dataset were docked on each of the cluster representatives of the ACE2-Spike. The phytochemicals identified belonged to Withania somnifera, Swertia chirayita, Tinospora cordifolia and Rutin DAB10, fulvestrant, elbasvir from FDA. Supplementary Information: The online version contains supplementary material available at 10.1007/s11696-021-01680-1.

Destabilization potential of beta sheet breaker peptides on Abeta fibril structure: an insight from molecular dynamics simulation study.

Alzheimer's disease is characterized by amyloid-ß aggregation. Currently, all the approved medications are to treat the symptoms but there is no clinically approved treatment for the cure or to prevent the progression of Alzheimer's disease (AD). Earlier reports suggest the use of small molecules and peptides to target and destabilize the amyloid fibril. The use of Beta Sheet Breaker (BSB) peptides seems to be a promising and attractive therapeutic approach as it can strongly bind and destabilize the preformed amyloid fibril. There are experimental studies describing the destabilization role of various BSB peptides, but the exact mechanism remains elusive. In the current work, an attempt is made to study the destabilization mechanism of different BSB peptides on preformed amyloid protofibril using molecular docking and simulations. Molecular docking of eight different BSB peptides of varying length (5-mer to 10-mer) on the Abeta protofibril was done. Docking was followed by multiple sets of molecular simulations for the Abeta protofibril-BSB peptide complex for each of the top ranked poses of the eight BSB peptides. As a control, multiple sets of simulations for the Abeta protofibril (APO) were also carried out. An increase in the RMSD, decrease in the number of interchain hydrogen bonds, destabilization of important salt bridge interactions (D23-K28), and destabilization of interchain hydrophobic interactions suggested the destabilization of Abeta protofibril by BSB peptides. The MM-GBSA free energy of binding for each of the BSB peptides was calculated to measure the binding affinity of BSB peptides to Abeta protofibril. Further residue wise contribution of free energy of binding was also calculated. The study showed that 7-mer peptides tend to bind strongly to Abeta protofibril as compared to other BSB peptides. The KKLVFFA peptide showed better destabilization potential as compared to the other BSB peptides. The details about the destabilization mechanism of BSB peptides will help in the design of other peptides for the therapeutic intervention for AD.

Drug repurposing studies targeting SARS-CoV-2: an ensemble docking approach on drug target 3C-like protease (3CLpro).

The COVID-19 pandemic has been responsible for several deaths worldwide. The causative agent behind this disease is the Severe Acute Respiratory Syndrome - novel Coronavirus 2 (SARS-CoV-2). SARS-CoV-2 belongs to the category of RNA viruses. The main protease, responsible for the cleavage of the viral polyprotein is considered as one of the hot targets for treating COVID-19. Earlier reports suggest the use of HIV anti-viral drugs for targeting the main protease of SARS-CoV, which caused SARS in the year 2002-2003. Hence, drug repurposing approach may prove to be useful in targeting the main protease of SARS-CoV-2. The high-resolution crystal structure of the main protease of SARS-CoV-2 (PDB ID: 6LU7) was used as the target. The Food and Drug Administration approved and SWEETLEAD database of drug molecules were screened. The apo form of the main protease was simulated for a cumulative of 150 ns and 10 µs open-source simulation data was used, to obtain conformations for ensemble docking. The representative structures for docking were selected using RMSD-based clustering and Markov State Modeling analysis. This ensemble docking approach for the main protease helped in exploring the conformational variation in the drug-binding site of the main protease leading to the efficient binding of more relevant drug molecules. The drugs obtained as top hits from the ensemble docking possessed anti-bacterial and anti-viral properties. This in silico ensemble docking approach would support the identification of potential candidates for repurposing against COVID-19.Communicated by Ramaswamy H. Sarma.

Remdesivir-bound and ligand-free simulations reveal the probable mechanism of inhibiting the RNA dependent RNA polymerase of severe acute respiratory syndrome coronavirus 2.

The efforts towards developing a potential drug against the current global pandemic, COVID-19, have increased in the past few months. Drug development strategies to target the RNA dependent RNA polymerase (RdRP) of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) are being tried worldwide. The gene encoding this protein, is known to be conserved amongst positive strand RNA viruses. This enables an avenue to repurpose the drugs designed against earlier reported inhibitors of RdRP. One such strong inhibitor is remdesivir which has been used against EBOLA infections. The binding of remdesivir to RdRP of SARS-CoV-2 has been studied using the classical molecular dynamics and ensemble docking approach. A comparative study of the simulations of RdRP in the apo and remdesivir-bound form revealed blocking of the template entry site in the presence of remdesivir. The conformation changes leading to this event were captured through principal component analysis. The conformational and thermodynamic parameters supported the experimental information available on the involvement of crucial arginine, serine and aspartate residues belonging to the conserved motifs in RdRP functioning. The catalytic site comprising of SER 759, ASP 760, and ASP 761 (SDD) was observed to form strong contacts with remdesivir. The significantly strong interactions of these residues with remdesivir may infer the latter's binding similar to the normal nucleotides thereby remaining unidentified by the exonuclease activity of RdRP. The ensemble docking of remdesivir too, comprehended the involvement of similar residues in interaction with the inhibitor. This information on crucial interactions between conserved residues of RdRP with remdesivir through in silico approaches may be useful in designing inhibitors.

TANGO: A high through-put conformation generation and semiempirical method-based optimization tool for ligand molecules.

Lead optimization is one of the crucial steps in the drug discovery pipeline. After identifying the lead molecule and obtaining its 2D geometry, understanding the best conformation it would attain in 3D still remains one of the most challenging steps in drug discovery. There have been multiple methods and algorithms that are directed toward achieving best conformation for the lead molecules. TANGO focuses on conformation generation and its optimization using semiempirical energy calculations. The conformation generation is based on torsion angle rotation of the exocyclic bonds. The energy calculations are performed using MOPAC. The unique feature of this tool lies in the implementation of Message Passing Interface (MPI) for conformation generation and semiempirical-based optimization. A well-defined architecture handling the input and output generation has been used. The master and slave approach to handle operations involved in torsion angle rotation and energy calculations has helped in load balancing the process of conformation generation. The benchmarking results suggest that TANGO scales significantly well across eight nodes with each node utilizing 16 cores. This tool may prove to very useful in high throughput generation of semiempirically optimized small molecule conformations. The use of semiempirical methods for optimization generates a conformational ensemble thereby helping to obtain stable and alternate stable conformers for a given ligand molecule. © 2018 Wiley Periodicals, Inc.

Detecting early stage structural changes in wild type, pathogenic and non-pathogenic prion variants using Markov state model.

The conversion of prion protein from normal to scrapie followed by the aggregation and deposition of this scrapie form leads to various neurodegenerative diseases. A few studies carried out by researchers suggest that E219K prion mutant (glutamate to lysine mutation at residue position 219) is more stable than wild type protein. However a similar point mutation E200K (glutamate to lysine mutation at residue position 200) is pathogenic. In this study we have carried out detailed atomistic simulation of the wild type, pathogenic mutant E200K and E219K mutant which provides more stability. The aim of the study was to detect the early structural changes present in all the three variants which might be responsible for the stability or for their conversion from PrPC to PrPSc. MSM based analyses have been carried out to find out the differences between WT, E200K and E219K systems. Markov state model (MSM) analysis was able to predict the intermediate states which helped to understand the effect of same mutation at two different locations. The MSM analysis was able to show that the extra stability of E219K mutant may be a result of the increase in number of native contacts, strong salt bridges and less random motions. While pathogenicity of E200K mutant can be attributed to loss of some crucial salt-bridge interactions, increased random motions between helix 2 and helix 3.

Correction to: Histone isoform H2A1H promotes attainment of distinct physiological states by altering chromatin dynamics.

After publication of this article [1], it was noticed Duane Smoot and Hassan Ashktorab who made and provided the cell line HFE145 were not included in the author list.

Acetylcholinesterase and Aß Aggregation Inhibition by Heterometallic Ruthenium(II)-Platinum(II) Polypyridyl Complexes.

Two heteronuclear ruthenium(II)-platinum(II) complexes [Ru(bpy)2(BPIMBp)PtCl2]2+ (3) and [Ru(phen)2(BPIMBp)PtCl2]2+ (4), where bpy = 2,2'-bipyridine, phen = 1,10-phenanthroline, and BPIMBp = 1,4'-bis[(2-pyridin-2-yl)-1H-imidazol-1-ylmethyl]-1,1'-biphenyl, have been designed and synthesized from their mononuclear precursors [Ru(bpy)2(BPIMBp)]2+ (1) and [Ru(phen)2(BPIMBp)]2+ (2) as multitarget molecules for Alzheimer's disease (AD). The inclusion of the cis-PtCl2 moiety facilitates the covalent interaction of Ru(II) polypyridyl complexes with amyloid ß (Aß) peptide. These multifunctional complexes act as inhibitors of acetylcholinesterase (AChE), Aß aggregation, and Cu-induced oxidative stress and protect neuronal cells against Aß-toxicity. The study highlights the design of metal based anti-Alzheimer's disease (AD) systems.

Biochemical and Biophysical Characterisation of Higher Oligomeric Structure of Rat Nucleosome Assembly Protein 1.

Nucleosome assembly protein 1 (NAP1) is a histone chaperone that exchanges histone H2A-H2B dimer from chromatin templates. Studies with yeast NAP1 (yNAP1) have revealed its existence as multiple oligomeric species in solution. Here, rat NAP1 (rNAP1), which is 98% identical to the human NAP1 (hNAP1) was used as a model to characterize the oligomeric structures of this protein in higher eukaryotes. Gel filtration chromatography and Dynamic light scattering of recombinant rNAP1 indicated that the protein exists as a complex mixture of multimeric species even at 500 mM ionic strength. The solution-state complexity remains unchanged even at higher ionic strengths. Equilibrium unfolding (ΔG 14.6 kcal mol- 1) shows that rNAP1, both dimeric and oligomeric, follow the two-state model of unfolding with no detectable intermediates. Homology modelling suggests that rat and yeast NAP1 share an overall similar structure with conserved domains. However, dissimilar substitutions like threonine and lysine with glycine in the ß-hairpin involved in oligomerization, possibly leads to the observed differences in the oligomerization propensity of the two proteins. Molecular dynamic simulation (MDS) of the two structures also revealed that rNAP1 dimer is more stable owing to the extensive hydrogen bonding in comparison to yNAP1. Further, in vitro kinase assay showed that the phosphorylation of rNAP1 favors oligomerization with no effect on its histone binding capacity. Our results clearly suggest that there are differences in the in-solution behavior of rNAP1 compared to yNAP1 which may have in vivo functional implications for the regulation of these complexes during chromatin assembly and rearrangement.

Histone isoform H2A1H promotes attainment of distinct physiological states by altering chromatin dynamics.

BACKGROUND: The distinct functional effects of the replication-dependent histone H2A isoforms have been demonstrated; however, the mechanistic basis of the non-redundancy remains unclear. Here, we have investigated the specific functional contribution of the histone H2A isoform H2A1H, which differs from another isoform H2A2A3 in the identity of only three amino acids. RESULTS: H2A1H exhibits varied expression levels in different normal tissues and human cancer cell lines (H2A1C in humans). It also promotes cell proliferation in a context-dependent manner when exogenously overexpressed. To uncover the molecular basis of the non-redundancy, equilibrium unfolding of recombinant H2A1H-H2B dimer was performed. We found that the M51L alteration at the H2A-H2B dimer interface decreases the temperature of melting of H2A1H-H2B by ~ 3 °C as compared to the H2A2A3-H2B dimer. This difference in the dimer stability is also reflected in the chromatin dynamics as H2A1H-containing nucleosomes are more stable owing to M51L and K99R substitutions. Molecular dynamic simulations suggest that these substitutions increase the number of hydrogen bonds and hydrophobic interactions of H2A1H, enabling it to form more stable nucleosomes. CONCLUSION: We show that the M51L and K99R substitutions, besides altering the stability of histone-histone and histone-DNA complexes, have the most prominent effect on cell proliferation, suggesting that the nucleosome stability is intimately linked with the physiological effects observed. Our work provides insights into the molecular basis of the non-redundancy of the histone H2A isoforms that are being increasingly reported to be functionally important in varied physiological contexts.

Ruthenium(II) polypyridyl complexes with hydrophobic ancillary ligand as Aß aggregation inhibitors.

The synthesis, spectral and electrochemical characterization of the complexes of the type [Ru(NN)2(txbg)](2+) where NN is 2,2'-bipyridine (bpy) (1), 1,10-phenanthroline (phen) (2), dipyrido [3,2-d:2',3f] quinoxaline (dpq) (3), and dipyrido[3,2-a:2',3'-c]phenazine (dppz) (4) which incorporate the tetra-xylene bipyridine glycoluril (txbg) as the ancillary ligand are described in detail. Crystal structures of ligand txbg and complex 2 were solved by single crystal X-ray diffraction. Thioflavin T (ThT) fluorescence and Transmission Electron Microscopy (TEM) results indicated that at micromolar concentration all complexes exhibit significant potential of Aß aggregation inhibition, while the ligand txbg displayed weak activity towards Aß aggregation. Complex 1 showed relatively low inhibition (70%) while complexes 2-4 inhibited nearly 100% Aß aggregation after 240 h of incubation. The similar potential of complexes 2-4 and absence of any trend in their activity with the planarity of polypyridyl ligands suggests there is no marked effect of planarity of coligands on their inhibitory potential. Further studies on acetylcholinesterase (AChE) inhibition indicated very weak activity of these complexes against AChE. Detailed interactions of Aß with both ligand and complex 2 have been studied by molecular modeling. Complex 2 showed interactions involving all three polypyridyl ligands with hydrophobic region of Aß. Furthermore, the toxicity of these complexes towards human neuroblastoma cells was evaluated by MTT assay and except complex 4, the complexes displayed very low toxicity.

Traversing the folding pathway of proteins using temperature-aided cascade molecular dynamics with conformation-dependent charges.

Protein folding is a multi-micro second time scale event and involves many conformational transitions. Crucial conformational transitions responsible for biological functions of biomolecules are difficult to capture using current state-of-the-art molecular dynamics (MD) simulations. Protein folding, being a stochastic process, witnesses these transitions as rare events. Many new methodologies have been proposed for observing these rare events. In this work, a temperature-aided cascade MD is proposed as a technique for studying the conformational transitions. Folding studies for Engrailed homeodomain and Immunoglobulin domain B of protein A have been carried out. Using this methodology, the unfolded structures with RMSD of 20 Å were folded to a structure with RMSD of 2 Å. Three sets of cascade MD runs were carried out using implicit solvation, explicit solvation, and charge updation scheme. In the charge updation scheme, charges based on the conformation obtained are calculated and are updated in the topology file. In all the simulations, the structure of 2 Å was reached within a few nanoseconds using these methods. Umbrella sampling has been performed using snapshots from the temperature-aided cascade MD simulation trajectory to build an entire conformational transition pathway. The advantage of the method is that the possible pathways for a particular reaction can be explored within a short duration of simulation time and the disadvantage is that the knowledge of the start and end state is required. The charge updation scheme adds the polarization effects in the force fields. This improves the electrostatic interaction among the atoms, which may help the protein to fold faster.